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Abstract 

We revise a recently proposed mechanism for spontaneous CP breaking at finite 
temperature in the Higgs sector of the Minimal Supersymmetric Standard Model, based 
on the contribution of squarks, charginos and neutralinos to the one- loop effective po- 
tential. We have included plasma effects for all bosons and added the contribution 
of neutral scalar and charged Higgses. While the former have little effect, the latter 
provides very strong extra constraints on the parameter space and change drastically 
the previous results. We find that CP can be spontaneously broken at the critical tem- 
perature of the electroweak phase transition without any fine-tuning in the parameter 
space. 
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1. The issue of spontaneous CP breaking (SCPB) at zero temperature has been recently 
proposed in supersymmetric theories as an appealing alternative |1| to the Peccei-Quinn 
mechanism to solve the strong CP problem, and also as a natural explanation for the small- 
ness of CP violating phases contributing to the neutron electron dipole moment (NEDM). 
In the minimal supersymmetric standard model (MSSM), SCPB at zero temperature can be 
triggered by radiative corrections ||. Unfortunately, as established on general grounds by 
Georgi and Pais || , it requires the existence of a Higgs boson with a mass of a few GeV || , 
which has been ruled out at LEP || . 

On the other hand, CP violation at finite temperature is one of the key ingredients to 
generate baryon asymmetry at the electroweak phase transition (7j. It has been recently 
realized || that temperature effects can trigger SCPB in the MSSM at the critical temper- 
ature of the electroweak phase transition, and baryon asymmetry with the help of tiny CP 
violating phases giving rise to a NEDM well below its present experimental limit M. The 
analysis of Refs. f|, fl was based on the contribution of squarks, charginos and neutralinos 
to the one-loop effective potential at finite temperature. 

In this paper we have analyzed SCPB in the MSSM at finite temperature, taking into 
account plasma effects for all bosons and including the contribution of the Higgs sector 
to the effective potential. We have found that plasma effects have little influence on the 
SCPB, but the contribution from neutral scalar and charged Higgses changes drastically 
the results found in Ref. || |9| giving rise to a much stronger constraint on the parameter 
space. Nevertheless we find that SCPB can be accomplished without any fine-tuning of the 
parameters. 



2. The most general scalar potential for the two-Higgs doublet model that is renormal- 
izable and SU(2) x [/(l) invariant is given by [TI| 

V = mx 2 !^! 2 + m 2 2 \H 2 \ 2 - {m z 2 E x E 2 + h.c.) + Ai|#i| 4 + A 2 |# 2 | 4 + A 3 |#i| 2 |# 2 | 2 

(1) 

+\ 4 \H 1 H 2 \ 2 + (A 5 (#i# 2 ) 2 + A 6 |#i| 2 #i# 2 + A 7 |# 2 | 2 #i# 2 + h.c), 

where m^ 2 , As, A6 and A7 are taken real, so that there is no explicit CP violation at tree level. 
However, SCPB can be triggered if the vacuum expectation values of the neutral components 
of the Higgs fields get a relative phase 6 7^ rnr, with integer n. (We can take (H°) = V\, 
(H%) = v 2 e tS without any loss of generality.) This requires, 

A 5 >0, 

-1 < cos s = !?SzMzM < L (2) 

4A 5 t>it; 2 

In the MSSM these conditions are not satisfied (at tree level) because in this case super- 
symmetry gives 
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where g and g' are the SU(2)l and U{1)y gauge couplings. Of course, after supersymmetry 
breaking, one—loop effects change these couplings and there are || finite contributions (pro- 
portional to the soft parameters) that can produce SCPB. At finite temperature, there will 
be additional corrections coming from the thermal modes running in the loops. 

We can write the one-loop effective potential, including the zero temperature terms, as 



b f 



(4) 



where g b (#/) denotes the number of degrees of freedom of bosons (fermions). In the bosonic 
part we sum over stops, neutral scalar and charged Higgses, and Goldstone bosons {b = 
tL,n, h, H, H ± ,G ± ), and in the fermionic one we consider neutralinos and charginos (/ = 
X ,^). Working in the 't Hooft-Landau gauge and in the DR renormalization scheme, we 
will decompose the bosonic components as 

V b = V nz {y 2 b ) + V-M)- (5) 
V nz includes the bosonic zero temperature part plus the finite temperature contribution 



coming from n^O Matsubara modes [11 
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where y\ = M. 2 /T 2 is the bosonic mass matrix rescaled with temperature, Q is the renor- 
malization scale, and the finite temperature one-loop term is given by the usual integral 
expresion with the n = part subtracted out, 
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The n = term is given by JTT 
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and now y 2 = (A4l+Ub)/T 2 , where lib is the thermal polarization matrix for bosons, contains 
the Debye masses flT2f . 

As the fermions do not have zero Matsubara modes, Vf is just given by fll| 



with 



V f = Tr 



647T 2 



log 



' T 2 y 2 j 

, Q 2 , 




dx x 2 log ( 1 + e 



-^Jx 2 +y 
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(10) 



and y 2 = A4 2 /T 2 , where Ai 2 is the fermionic mass matrix. 

The bosonic and fermionic mass matrices are functions of the classical fields H°, H^-, 
H\ and i?2°, and one can expand the potential (4) in powers of them to obtain the one- 
loop corrections to the tree level couplings of Eq. (1). Assuming as in Ref. |y equal soft 
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supersymmetry breaking masses for left and right-handed squarks 1 ttlq 
can easily find the contributions from stop loops to m|, A5, A6 and A7: 



raft = m one 

(11) 
(12) 

(13) 
(14) 



where Q fi is the supersymmetric Higgs mixing mass, A t the trilinear soft breaking parameter 
corresponding to h t Q • H 2 U in the superpotential, 
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(and a similar definition for higher derivatives), and now y~ = m 2 /T 2 , y~ = (m 2 + Uf)/T 2 
with 0, g 



(17) 



where g s is the SU(3) gauge coupling. We are neglecting in Eq. (17) contributions of order 
g 2 ,g >2 and h 2 compared to the g 2 part. In this way the field independent mass matrix for 
stops with the thermal screening included is proportional to the identity matrix and the 
expressions (11-14) apply. 

The mass matrices of neutral and charged Higgses at the origin, H° = iff = 0, are never 
proportional to the identity and so the corresponding one-loop corrections involve finite 
differences instead of derivatives of the If, functions. In particular, 
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1 This is just a technical simplification such that the mass matrix at the origin, = = 0, is propor- 
tional to the identity matrix. 

2 The high-T expansion of Eqs. (11-14) coincides with Eqs. (15-18) in Ref. ||, except that we find 1/512 
instead of 1/256 in Eq. (16) of [§. 
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where the contribution to Aw proportional to (g 2 + g' 2 ) inside the curly brackets comes 
from the neutral scalar Higgs sector, and the rest, proportional to g 2 , from the charged 
Higgs sector. Alh and A 3 Ih in Eqs. (18-21) are defined by 

lllh = o o r 
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and y\ = m^/T 2 , y± = m\jT 2 as usual, with m\ the Higgs masses for zero classical fields: 



m i 



- (ml + ml ± yj (m 2 , — m 2 ) 2 + 4m 2 j 



and m± the corresponding masses in the presence of screening: 



m i 



- (ml + ml ± ^/ (ml — m 2 ) 2 + 4m 2 j 



(23) 



(24) 



where m 2 and m\ are functions of m 2 , tan/? and the other parameters of the theory, by 
means of the zero-temperature (including the one-loop corrections) minimization conditions, 



and IO, IT1 
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mi + U hl =mi + ^(3g 2 + g' 2 )T 2 
ml + U h2 = m 2 2 + l(3g 2 + g' 2 + 6h 2 )T 2 . 



(25) 
(26) 



We can also assume that 1 /i 2 = M\ = Mf (where M2 and Mi are the soft Majorana 
masses for SU{2) L and U(l)y gauginos, respectively) and then find, for charginos: 
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3 We correct here Eq. (26) of Ref. || where the proportionality coefficient of A^A5,6,7 was written as 
(.g 2 +.g' 2 )/2 5 2 . 
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where now y"^ = f^/T 2 and all the integrals are fermionic. Derivatives are defined, as usual, 

h = v'f = ^ f v M)- ( 32 ) 



3. We analyze now the effect described in the previous section. The parameter m\ will be 
traded by m\, the physical mass of the pseudoscalar which includes all radiative corrections 
(at zero-temperature). None of conditions (2) is satisfied at tree level (in particular, A5 = 
and 5 = wr), neither after including one-loop radiative corrections at zero-temperature, if 
we fix rriA beyond its present experimental limit |J. However, they can be satisfied when 
radiative corrections at finite-temperature are included. In particular, A a , (a = 5, 6, 7) 
acquires a renormalized value at finite temperature 

A a (T)= AWA - ( 33 ) 

i=s,h,c,n 

as given by Eqs. (12-14), (19-21), (28-29) and (31). Also ml gets renormalized as 

ml(T)=ml+ J2 (34) 

i=s,h,c,n 

from Eqs. (11), (18), (27) and (30). We have computed, in (33) and (34), A^' c ' n ^ numerically, 
from the integrals (7) and (10), while A^ is computed using its finite temperature expansion, 
because in the Higgs sector the corresponding value of the integral (7) is well defined only 
through its high-T expansion. 

A glance at the previous section (Eqs. (12), (19), (28) and (31)) shows that A^As < 0, 
while A ( c > n >\ 5 > 0, so that the contribution of fermions (charginos and neutralinos) have to 
compensate that of squarks, if no other contributions were considered, so that one would 
obtain in the plane (/i, m) a lower bound contour, as in Ref . || . We have analyzed the 
effects of the Debye screening on the latter contour and found it gets shifted towards smaller 
(~ 5%) values of m. We can conclude from this that plasma corrections in the squark sector 
are not very important. Now the contribution from the Higgs sector (where plasma effects 
are essential to make the effective potential at H° = H 1 ^ = real) is negative, and has an 
infrared singularity (oc — (T/m_)), responsible for the failure of perturbative expansion |11 



for values of (/i, m) such that m_ = (long-dashed lines in Figs. 1-3). This creates in the 
plane (/z, m) an upper bound contour which severely constrains the region of the parameter 
space where A5 > 0. This behaviour is exemplified in Figs. 1-3 where A5 = is plotted in the 
plane (/1, m) for different values of T, rriA and the supersymmetric parameters 0. We can 
see from Figs, lb, 2b and 3b that the lower and upper contours provide a combined region 
where A5 > 0. In the absence of mixing A t = 0, Figs, la, 2a and 3a, there is no contribution 
coming from the squark sector and, consequently, the lower contour disappears. Only the 
upper contour (from the Higgs sector) constrains the region where A5 > 0. If we compare 

4 We have checked that the theory remains perturbative, i.e. that (3± = g 2 T/fh± < 1, in the region where 
A 5 > in Figs. 1-3. 
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the region where As(T) > for different temperatures, e.g. Figs. 1 and 2, we see it decreases 
with decreasing T, and for any given point in the plane (/i, in) the condition A5 > will stop 
being satisfied at a given temperature. 

In the region where As(T) > the condition | cos<5(T) |< 1 where 

cos 6(T) - ml(T) - X 6 (T)vl(T) - X 7 (T)v!(T) 

is also required to have SCPB. This condition is satisfied for values of v (T) (v 2 (T) = v\ (T) + 
f|(T)) such that V-(T) < v(T) < v + (T), where 

v±(T) = ^jjL {A 6 (T) + A 7 (T) tan 2 (3(T) T 4A 5 (T) tan/3(T)} _1 (36) 

For a given point in the plane (/i, in) the values of v±(T) are very close to each other and 
depend on tan/3(T). For a first-order phase transition the values of tan/3(T) and (v(T)) at 
the critical temperature are dynamically fixed by all the parameters of the supersymmetric 



theory. In the absence of a complete analysis of the phase transition in the MSSM JLJ], we 
will take tan/3(T) = tan/3, which we know is a good approximation at least for large values 
of rriA, to get an estimate of the values of v± for the different points in the plane (//, in). At 
very high temperatures, (v(T)) = and the condition (2) will not be satisfied in general for a 
fixed point in the plane (//, in). At the critical temperature T c the field will go from the origin 
to (v(T c )) and produce a maximal CP violation (AS = ±7r) if < f_(T c ) < v + (T c ) < (v(T c )) 
for the chosen point in the plane in). 

We have plotted in Figs. 1-3 level contours corresponding to different values of < 
V -(T) < T. If the phase transition is strong enough, i.e. (v(T)) ~ T, then the whole region 
with V-.(T) < T and A5 > satisfies the condition of CP violation. If the phase transition 
is stronger (weaker), i.e. (v(T)) > T ((v(T)) < T), then the allowed region showed in Figs. 
1-3 will be further enlarged (reduced). For values of (tan/3(T)) different from tan/5, the 
allowed regions in Figs. 1-3 would get distorted, but still a wide region would appear. 



4. In conclusion we have estimated the region of maximal CP violation at the critical 
temperature of the electroweak phase transition in the space of supersymmetric parameters. 
We have studied for simplicity the degenerate case of uiq = mu in the stop sector, and 
Mi = M 2 = n in the gaugino/higgsino sector. We have found that the Higgs sector changes 
dramatically previous results on the same effect. Nevertheless the mechanism can work in 
a wide region of the parameter space without any fine-tuning, and trigger easily maximal 
CP violation during the first order phase transition. However this requires, as a general 
tendency, small values of and large values of tan/3, while preliminary results on the 
phase transition in the MSSM [|13|, [14| seem to point towards the opposite direction in order 
not to wash out the generated baryon asymmetry. The analysis of the non-degenerate case 
itlq ^ mu, Mi ^ M2 7^ yU (now in progress) seems necessary to find out whether this 
mechanism can be realistically used to generate baryon asymmetry at the phase transition 
in the MSSM. 
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Fig.l la) Contours of m_ = (long-dashed line), A5 = (horizontal solid line), v-(T)/T = 
1 (outer vertical solid lines) and v-(T)/T = (inner vertical solid lines) for m t = 
160 GeV, T = 150 GeV, m A = 50 GeV, tan/5 = 8 and A t = 0; lb) The same as in la) 
but for A t = 50 GeV. Here the contour of A5 = is closed on the left and the contours 
of v-(T)/T = 1 (outer vertical solid lines) and v-(T)/T = 0.6 (short-dashed line) are 
shown. 

Fig.2 The same as in Fig. 1 but with T = 100 GeV and tan (3 = 18. In 2b) the short-dashed 
line corresponds to v_(T)/T = 0.7 

Fig. 3 The same as in Fig. 1 but for m A = 75 GeV and tan (3 = 18. In 3b) the short-dashed 
line corresponds to v-(T)/T = 0.5 
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